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Abstract 

Water  management  is  critically  important  for  polymer  electrolyte  membrane  fuel  cells  (PEMFC),  and  is  complicated  by  the  electroosmotic 
flow  of  water  from  anode  to  cathode  through  the  polymer  electrolyte  membrane.  In  this  study,  electroosmotic  flow  in  polymer  electrolyte 
membranes  is  modeled  incorporating  the  electrokinetic  effect,  and  key  parameters  affecting  the  PEM  fuel  cell  performance  are  identified. 
The  governing  Poisson-Boltzmann  and  the  Navier-Stokes  equations  were  solved  numerically  for  a  single  membrane  pore  to  determine  the 
electroosmotic  flow  through  the  membrane  over  a  wide  range  of  geometrical  and  operating  conditions.  It  was  found  that  the  electroosmotic 
drag  coefficient,  ^fDrag,  increases  with  the  pore  diameter.  The  membrane  thickness  has  a  significant  effect  on  the  electroosmotic  flow.  At 
constant  cell  voltage,  the  electroosmotic  flow  through  thicker  membranes  (e.g.  Nafion  117)  is  reduced  because  of  the  reduced  electric  field 
strength.  The  pressure  difference  required  to  stop  electroosmotic  flow  is  very  large  due  to  the  extremely  small  pore  diameters.  In  the  presence 
of  sulphuric  acid,  numerical  results  have  revealed  that  the  electroosmotic  flow  increases  with  the  acid  concentration. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Polymer  electrolyte  membrane  fuel  cell  (PEMFC)  is 
one  of  the  most  promising  candidates  as  zero  emission 
power  sources  for  both  stationary  and  mobile  applications. 
It  is  highly  efficient  and  environmentally  clean.  During  the 
PEMFC  operations,  the  polymer  electrolyte  membrane  pro¬ 
vides  the  transport  medium  for  the  protons,  generated  during 
the  anodic  oxidation  of  the  fuel,  to  migrate  to  the  cathode 
and  complete  the  electric  and  mass  transport  circuits.  There¬ 
fore,  high  proton  conductivity  in  the  membrane  is  essential 
to  achieve  high  performance  in  PEM  fuel  cells.  The  protonic 
current  through  the  membrane,  however,  tends  to  drag  wa¬ 
ter  molecules  along  resulting  in  the  water  depletion  at  the 
anode  side  and  water  accumulation  (hence  flooding)  at  the 
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cathode.  This  phenomenon,  also  called  electroosmotic  flow 
(or  drag),  increases  the  membrane  resistance  and  reduces  fuel 
cell  performance  [1,2]. 

The  water  content  in  the  solid  polymer  electrolyte  mem¬ 
brane  is  determined  by  the  balance  of  water  production  at 
the  cathode  and  its  transport  during  the  fuel  cell  operation. 
Fig.  1  depicts  schematically  the  mechanisms  of  water  trans¬ 
port  in  a  PEM  fuel  cell.  Water  is  carried  into  the  fuel  cell 
via  the  humidified  gas  streams  entering  gas  diffusion  elec¬ 
trodes.  Some  combination  of  water  vapor  and  liquid  water 
passes  through  each  electrode  to  the  electrode/membrane  in¬ 
terface.  Water  crossing  this  interface  assists  in  the  hydration 
of  the  membrane.  An  additional  source  of  water  is  the  oxy¬ 
gen  reduction  reaction  occurring  at  the  cathode.  Inside  the 
membrane,  water  is  transported  in  two  ways:  electroosmotic 
flow  and  back  diffusion  from  cathode  to  anode.  While  both 
electroosmotic  drag  and  water  production  due  to  oxygen  re¬ 
duction  tend  to  create  an  excess  of  water  near  the  cathode, 
back  diffusion  of  water  from  cathode  to  anode  tends  to  op¬ 
pose  the  build-up  of  water  and  flatten  the  water  concentration 
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Nomenclature 

a 

pore  radius  (m) 

e 

fundamental  electronic  charge  (1.602  x 

i<r19c) 

E 

electric  field  strength  (V  m-1) 

£* 

non-dimensional  term,  Eq.  (14) 

n 

non-dimensional  electroosmotic  force,  Eq. 
(13) 

Fp 

non-dimensional  pressure  force,  Eq.  (13) 

i 

ion  number  in  the  electrolyte 

Kn 

Knudsen  number 

&B 

Boltzmann  constant  (1.381  x  10-23  JK-1) 

^Drag 

electroosmotic  drag  coefficient 

L 

pore  length  (m) 

M 

molecular  weight  (kg  mol-1) 

M* 

non-dimensional  term,  Eq.  (14) 

n 

iteration  number 

no 0 

bulk  ionic  number  concentration  (m-3) 

N 

local  number  concentration  (m-3) 

N 

molar  flow  rate  (molm  -3  s  -1) 

Na 

Avogadro  number  (6.02  x  1023  parti¬ 
cles  mol-1) 

P 

pressure  (Pa) 

Q 

volumetric  flow  rate  (m3  s-1) 

r 

pore  radius  (m) 

r* 

non-dimensional  pore  radius 

T 

temperature  (K) 

^EO 

electroosmotic  velocity  (ms-1) 

u 

reference  velocity  (e.g.  1  mm  s-1) 

u* 

non-dimensional  velocity 

* 

us 

non-dimensional  slip  velocity 

W 

wall 

z 

valance  of  the  dissolved  ions 

Greek 

letters 

€ 

dimensionless  dielectric  constant 

eo 

permittivity  of  vacuum  (8.85  x  10-12 

C2  N-1m-2) 

K 

Debye-Huckle  parameter  (m-1) 

P 

viscosity  (Pa  s) 

\jf 

potential  (V) 

non-dimensional  potential 

p 

density  (kgm-3) 

Pe 

net  volumetric  charge  density  (C  m-3) 

5 

lattice  spacing  (m) 

A 

difference 

T] 

distance  from  the  pore  wall  (A) 

zeta  potential  (V) 

profile  across  the  membrane.  It  is  therefore  very  important  to 
determine  and  manage  the  extent  of  the  electroosmotic  flow 
in  the  membrane  to  optimize  the  PEM  fuel  cell  performance. 


ANODE  CATHODE 

H2  (g)  —  2H++  2e  02(g)  +  4H*+4e  —  2H2P 


Fig.  1.  Schematic  of  water  transport  and  electroosmotic  flow  through  the 
polymer  electrolyte  membrane  in  a  PEM  fuel  cell. 

The  electroosmotic  drag  coefficient,  Kvmg,  defined  as 
the  number  of  water  molecules  transferred  through  the 
membrane  per  proton  in  the  absence  of  other  water  transport 
mechanisms,  has  been  measured  by  a  variety  of  techniques. 
The  measured  value  of  this  coefficient  for  Nation,  the  most 
commonly  used  membrane,  however,  exhibits  a  wide  scatter. 
LaConti  et  al.  [3]  measured  the  transferred  water  volume  us¬ 
ing  capillary  tubes.  They  reported  a  linear  increase  in  Kj^mg 
from  0  (in  the  dry  state)  to  4-5  (in  the  fully  hydrated  state)  for 
Nation  membranes  with  equivalent  weights  of  1150-1275 
and  found  no  temperature  dependence.  Zawodzinski  et  al.  [1] 
used  a  similar  method  and  obtained  Kj)mg  between  0.9  and 
2.5  for  partially  and  fully  hydrated  Nation  117  membrane 
at  300  K,  respectively.  Fuller  and  Newman  [4]  measured 
electroosmotic  drag  for  Nation  117  in  a  H2O  concentration 
cell.  They  reported  K^mg  ^  1.4  for  hydrated  membrane 
and  a  decrease  to  Kj)mg  =  0  for  the  dry  membrane  at  25  °C. 
From  the  same  kind  of  measurement,  Zawodzinski  et  al.  [2] 
later  found  AT>rag  =  1  for  vapor-equilibrated  membranes 
over  a  wide  range  of  water  content  and  almost  independent 
of  the  membrane  micro  structure.  In  fact,  this  coefficient 
would  best  represent  the  macroscopic  phenomenon  of 
electroosmotic  drag  in  an  operating  PEM  fuel  cell,  since 
the  majority  of  the  water  molecules  in  the  membrane  are  in 
the  form  of  vapor  at  the  fuel  cell  operating  conditions.  Xie 
and  Okada  [5]  measured  K& rag  =  2.6  for  a  fully  hydrated 
membrane  using  the  streaming  potential.  Recently,  Ren  et  al. 
[6]  presented  a  method  of  measuring  electroosmotic  drag  in 
a  direct  methanol  fuel  cell  and  obtained,  for  a  fully  hydrated 
membrane,  an  increase  of  K& rag  as  a  function  of  temperature 
from  1.9  at  15  °C  to  5.1  at  130  °C.  The  wide  range  of  scatter 
in  the  reported  K& rag  clearly  shows  the  sensitivity  and 
difficulties  involved  with  the  experimental  techniques. 

The  PEM  electroosmotic  driven  flow  in  general,  and  the 
electroosmotic  drag  in  particular,  can  also  be  investigated 
with  the  knowledge  of  the  pore  geometry  and  electrokinetic 
effect.  It  is  therefore  the  aim  of  this  study  to  develop  a 
model  that  describes,  under  steady- state  conditions,  the 
electroosmotic  flow  in  an  operating  fuel  cell  using  the 
fundamental  physical  mechanisms.  The  formulation  is  from 
a  macroscopic  point  of  view,  but  with  sufficient  incorporation 
of  the  nanoscopic  information,  although  it  may  be  much 
simplified,  such  that  the  resulting  model  has  sufficient  pre- 
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dictive  capability  with  acceptable  accuracy  for  macroscopic 
engineering  application.  Therefore,  it  is  not  interested  in 
including  all  the  fine  features  of  the  nanoscopic  or  molecular 
level  mechanisms  of  the  proton  transport  processes.  To 
this  end,  the  Poisson-Boltzmann  and  the  Navier-Stokes 
equations  were  solved  numerically  to  investigate  the  elec- 
troosmotic  flow  in  a  simplified  PEM  pore  configuration  over 
a  wide  range  of  geometrical  and  operating  conditions. 


2.  Model  formulation 

2.7.  The  physical  problem 

Perfluorinated  sulphonic  membranes,  such  as  Nation,  are 
the  ionomeric  materials  typically  used  in  PEM  fuel  cells. 
The  study  of  water  and  ions  transport  in  polymer  electrolyte 
membranes  requires  extensive  knowledge  of  the  membrane 
morphology,  water  sorption  thermodynamics  and  proton  dis¬ 
tribution  in  the  membrane,  interactions  between  the  protons 
and  the  membrane,  and  proton  transport  in  aqueous  solution. 
A  great  number  of  experimental  techniques  have  been  applied 
in  an  effort  to  understand  the  morphology  and  ion  distribu¬ 
tion  of  polymer  electrolyte  membranes  and  the  correlation  of 
these  two  properties.  Some  of  these  techniques  include  small 
and  wide  angle  X-ray  scattering  [7,9],  dielectric  relaxation 
[10],  dielectric  spectroscopy  [11],  electron  microscopy  [12], 
and  AC  impedance  investigations  [1 3] ,  to  quote  a  few.  Despite 
the  substantial  work  conducted  for  characterization  of  poly¬ 
mer  electrolyte  membranes,  a  fundamental  understanding  of 
the  mechanisms  of  proton  and  water  transport  is  not  fully 
understood  [14].  This  is  largely  due  to  the  very  inhomoge¬ 
neous  nature  of  these  membranes  when  hydrated.  However, 
majority  of  these  studies  consider  polymer  electrolyte  mem¬ 
brane  as  being  made  of  three  microphase- separated  regions 


A 


t~ Flurocarbon  backbone  material 
0  -S03  0  SOt  ©  H*  O  H20 


Fig.  2.  Nafion  cluster  network  model  [7,8]. 


as  illustrated  in  Fig.  2.  According  to  these  models,  region  A 
consists  of  hydrophobic  fluorocarbon  backbone  material.  The 
ion-cluster  regions,  C,  generally  described  as  pores  or  aque¬ 
ous  solution  environments,  with  sulphonate  groups  (-SO3-) 
as  fixed-charge  sites  and  movable  ions  (usually  H+  for  acidic 
membranes)  as  counterions.  The  interconnecting  region  B 
is  an  amorphous  hydrophobic  region  of  lower  ionic  content. 
The  hydrophilic  ion  clusters  C  and  region  B  are  responsible 
for  ionic  conduction.  A  network  of  these  hydrophilic  clusters 
of  the  size  50-60  A  is  cross-connected  with  other  hydrophilic 
channels  of  smaller  diameters  10-20  A  [8]. 

A  number  of  physical  models  have  also  been  developed 
to  identify  the  mechanisms  of  water  and  ion  transport  in 
the  ionomer  membranes.  The  hydrodynamic  model  approach 
[15],  solubility  and  transport  models  [16],  partition  coeffi¬ 
cient  model  [17,18],  and  random  network  model  of  a  micro- 
porous  polymer  electrolyte  membrane  [19]  are  only  a  few 
examples.  These  models  do  not  represent  the  true  morpho¬ 
logical  features  associated  with  the  ionomer  membranes,  and 
rely  on  the  input  parameters.  As  a  result,  their  predictions 
often  deviate  from  the  experimental  measurements. 

Molecular  theoretical  studies  [20,21]  directed  towards 
understanding  these  remarkable  properties  of  ionomer 
membranes  are  limited  for  several  reasons.  Polar-solvent 
environments  profoundly  affect  charge  transport  processes 
like  proton  transfer,  where  collective  solvent  effects,  in¬ 
cluding  orientational  polarization,  require  consideration 
of  large  molecular  ensembles  in  computer  simulations. 
Proton  motion  exhibits  quantum  mechanical  behavior  such 
as  zero-point  energy  constraints  and  tunneling,  requiring 
quantum  mechanical  treatment  that  is  computationally 
challenging  and  expensive  in  multi-body  systems. 

On  the  other  hand,  the  proton  transport  through  the 
membrane  electrolyte  has  been  conventionally  modeled  by 
the  phenomenological  Ohm’s  law  in  most  of  engineering- 
oriented  calculations  for  operating  fuel  cells  [22].  This 
macro-homogeneous  approach  does  not  offer  any  fundamen¬ 
tal  insight  on  the  mechanism  of  actual  proton  transport  pro¬ 
cesses  at  the  molecular  level.  The  present  model  envisions 
the  electroosmotic  drag  phenomenon  from  another  perspec¬ 
tive.  The  hydrophilic  tortuous  electrolyte  regions  within  the 
membrane  are  considered  as  continuous,  parallel  conduction 
pathways  with  average  diameters  ranging  from  20  to  80  A, 
perpendicular  to  the  membrane  surface  with  uniform  sul¬ 
phonic  ions  fixed  at  the  wall  resulting  in  an  equipotential 
surface.  The  physical  and  electrochemical  properties  (with 
the  exception  of  electrolyte  permittivity)  are  assumed  to  be 
constant  and  independent  of  the  position.  This  approach  is 
advantageous  because  it  is  physically  sound,  simple,  and  re¬ 
quires  no  experimental  data.  As  a  result,  it  can  be  used  for 
parametric  study  and  provides  some  design  and  optimization 
rules  without  going  into  the  detailed  molecular  level  transport 
phenomena. 

The  presence  of  the  negative  sulphonic  charge  on  the  wall 
results  in  the  formation  of  a  diffuse  layer  of  counterions  (e.g. 
protons  in  the  present  case)  where  the  concentration  near  the 
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ANODE 


CATHODE 


©  Sulfonic  group  6-  Protons  C  Water  molecules 


Electric  Field 


Fig.  3.  Electroosmotic  flow  of  water  in  a  membrane  pore  of  a  PEM  fuel  cell. 

sulphonic  wall  is  at  the  highest  and  decreases  gradually  with 
distance  until  it  reaches  equilibrium  [23].  The  resulting  elec¬ 
tric  potential  is  therefore  minimum  at  the  wall  and  gradually 
increases  towards  an  equilibrium  level  as  the  concentration 
of  protons  approaches  the  bulk  concentration  at  farther  dis¬ 
tances. 

If  an  electric  field  is  applied  across  the  membrane,  as  in  an 
operating  fuel  cell,  an  electrokinetic  body  force  is  exerted  on 
the  charged  analytes  (e.g.  free  H+  ions  not  interacting  with 
S03").  The  protons  will  then  move  in  the  narrow  hydrophilic 
pores  under  the  influence  of  the  applied  electric  field,  taking 
with  them  the  hydrating  water  molecules,  and  resulting  in  an 
electroosmotic  flow  as  illustrated  in  Fig.  3.  The  fluid  move¬ 
ment  is  carried  through  to  the  rest  of  the  fluid  in  the  pores  by 
viscous  forces.  It  should  be  pointed  out  that  the  electric  po¬ 
tential  (or  field)  across  the  membrane,  as  shown  in  Fig.  3,  is 
after  considering  the  effect  of  the  electric  double  layers  near 
the  electrode  and  the  membrane  electrolyte  interface  [24] . 

2.2.  The  charge  distributions 

The  electric  potential  distribution  existing  in  the  charged 
analytes  flowing  inside  a  membrane  pore  shown  in  Fig.  3  is 
due  to  the  presence  of  diffuse  layer  in  a  cylindrical  pore;  it 
can  be  related  to  the  local  volumetric  net  charge  density  of 
the  electrolyte,  pe,  using  theory  of  the  electrostatics  [25]: 

V  •  (eVi/r)  =  — —  (1) 

eo 


imated  mathematically  as: 


6  = 


'17 

<  17  +  23.707  -7.8) 
,  81 


rj  <  7.8  A 
7.8  <  o  <  10.5  A 
0  >  10.5  A 


(2) 


The  net  charge  density  can  be  expressed  in  terms  of  the  Boltz¬ 
mann  distribution  as  [25]: 


Pe  =  ^2 ZienUoo  exp  wj  (3) 

where  n  is  the  valence  of  the  dissolved  ions,  n ;?00  is  the  ionic 
number  concentration  in  the  bulk  solution,  e  is  the  fundamen¬ 
tal  electronic  charge,  k#  is  Boltzmann’s  constant,  and  T is  the 
absolute  temperature. 

Eq.  (3)  can  be  simplified  for  the  proton-conducting  mem¬ 
brane  used  in  PEM  fuel  cells  where  only  counterions  (H+, 
z  =  +1)  are  present: 


Pe  =  zen  oc  exp 


(4) 


Substitution  of  Eq.  (4)into  Eq.  (1)  and  transforming  it  into 
non-dimensional  form  results  in  the  Poisson-Boltzmann 
equation  as  follows: 


j_j_  /  ,dr\ 

r*  dr*  \  dr*  ) 


1 

2 


(k\o)2  exp (— ^*) 


(5) 


The  non-dimensional  variables  are  defined  as: 


\lr*  = 


ifr 

(kBT)/(ze) 


and  r* 


r 

a 


(6) 


where  r  is  the  radial  distance  and  a  is  the  pore  radius. 

The  so-called  Debye-Huckle  parameter,  k\,  which  is  a 
measure  of  the  diffuse  layer  thickness,  is  defined  as  [27]: 


/  2nooZ2e2\  1/2 

V  66o k^T  ) 


(7) 


Eq.  (5)  is  subject  to  the  following  boundary  conditions: 
diA* 


atr*  =0, 


dr* 


=  0;  and 


r*  =  1,  = 


tAw 


(kBT)/(ze) 


=  ^w 


(8) 


where  6  is  the  dimensionless  dielectric  constant  of  the  elec¬ 
trolyte,  is  the  electrostatic  potential,  and  6o  is  the  permit¬ 
tivity  of  vacuum.  The  dielectric  constant  of  the  electrolyte 
varies  as  a  function  of  the  radial  distance  from  the  pore  wall 
due  to  the  change  in  the  orientations  of  the  water  molecules. 
A  dielectric  saturation  model  was  developed  by  Paddison  and 
Paul  [26]  for  fully  hydrated  Nafion  membranes.  According 
to  this  model,  dielectric  constant  remains  that  of  bulk  water 
for  distances  of  at  least  12  A  from  the  wall,  but  begins  to  fall 
sharply  as  the  arrays  of  -SO3-  fixed  sites  are  approached. 
Based  on  their  model,  variations  of  dielectric  constant  as  a 
function  of  the  distance  from  the  pore  wall,  77,  can  be  approx- 


The  potential  at  the  pore  wall,  iAw>  which  is  also  the  bound¬ 
ary  of  the  fluid  flow  problem,  is  called  zeta  potential,  £.  This 
potential  is  typically  determined  empirically  from  electroos¬ 
motic  or  streaming  potential  flow  measurements  [27] . 

23.  The  velocity  distribution 

The  electroosmotic-induced  flow  can  be  obtained  with  the 
knowledge  of  the  liquid  velocity  distribution.  With  the  pore 
sizes  in  the  order  of  nanometers,  the  Reynolds  number  is 
very  small  and  the  entry  length  effects  can  be  neglected  such 
that  the  axial  velocity  would  be  a  function  of  radius  (r)  only. 
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However,  the  minute  pore  geometry  may  not  hold  enough 
number  of  molecules  to  satisfy  completely  the  concept  of 
continuum.  The  corresponding  Knudsen  {Kn)  number  can  be 
calculated  from: 

5 

Kn  =  -  (9) 

a 

where  the  lattice  spacing,  8 ,  is  about  3  A  for  water  [28].  For 
an  average  pore  size  of  50  A,  the  Knudsen  number  is  about 
0.06  and  a  slip  flow  regime  is  expected  to  exist  in  the  pore.  In 
the  absence  of  thermal  creep,  the  dimensionless  slip  velocity 
at  the  pore  wall  is  calculated  as: 

rk/* 

ul  =  Kn-±  (10) 

on 

where  n  and  subscript  s  stand  for  the  outward  normal  (unit) 
vector  and  the  tangential  (unit)  vector  [29] . 

For  the  pore  wall  negatively  charged,  the  electroosmotic 
flow  will  be  in  the  direction  towards  the  cathode  as  illus¬ 
trated  in  Fig.  3.  The  body  force  per  unit  volume  created 
by  electric  field  acting  on  the  charged  particles  is  given  by 
Fe  =  PqE,  where  E  is  the  electric  field  strength.  The  Navier- 
Stokes  equations  for  a  Newtonian  laminar  liquid  flow  with 
constant  density,  p,  and  viscosity,  /x,  in  a  horizontal  cylindri¬ 
cal  capillary  under  the  influence  of  both  an  external  electric 
force  and  pressure  gradient  can  be  expressed  as  [25]: 

/xV2w  —  Vp  +  PqE  =  0  (11) 


representing  the  balance  among  the  net  viscous  force,  pres¬ 
sure  force,  and  electric  force  acting  on  the  fluid  particles  for 
the  steady  flow  under  consideration. 

Eq.  (11)  can  be  non-dimensionalized  as: 


1  d 
r*  dr* 


+  F*  +  F*  =  0 


(12) 


where  w*  is  the  liquid  velocity  normalized  with  some  ref¬ 
erence  velocity,  U,  considered  as  lmms-1  throughout  this 
work.  Fp  and  are  non-dimensional  pressure  and  elec¬ 
troosmotic  body  forces,  respectively,  defined  as: 

A  Pa2 

F*P  =  ~,  Fe  =  M*  p*E*  (13) 

and 

M*  =  zencc^f  e*  =  — 

fiUL  ’  f 

p*  =  -15—  =  exp  (14) 

zen  oo 

where  A  P  represents  the  pressure  difference  between  the  an¬ 
ode  and  cathode  ends  of  the  pore  and  L  is  the  capillary  length. 
The  dimensionless  group,  M*,  is  the  ratio  of  the  electric  to 
viscous  forces  per  unit  volume.  Eq.  (12)  shows  that  pore  hy¬ 
drodynamics  is  influenced  by  two  driving  forces:  pressure 
difference  across  the  membrane  and  the  electric  force  acting 
on  the  charge  carriers. 


Eq.  (12)  is  subject  to  the  following  boundary  conditions: 

atr*  =  0,  =  0;  and 

dr* 

r*  =  1,  u*  =  ul  (slip  velocity)  (15) 


3.  Solution  method 

There  is  no  exact  solution  for  the  non-linear  Poisson- 
Boltzmann  equation,  Eq.  (5).  A  finite  volume  scheme  with 
variable  cell  spacing  was  developed  to  discretize  the  equa¬ 
tion.  To  facilitate  convergence,  the  exponential  term  has  been 
linearized  as  follows: 

exp(— fl+l)  =  exp (— i/r*)  -  exp(— VO  0*+i  -  ft)  (16) 

where  n  refers  to  the  iteration  number. 

The  resulting  system  of  algebraic  equations  was  solved  us¬ 
ing  the  Gauss-Seidel  iterative  technique.  Totally,  500  control 
volumes  were  created.  The  grid  spacings  were  reduced  suc¬ 
cessively  by  5%  to  condense  the  control  volume  sizes  in  the 
region  near  the  wall,  where  the  electroosmotic  effect  is  sig¬ 
nificant.  The  calculated  electric  potential  distribution  is  then 
substituted  into  the  momentum  equation,  Eq.  (12),  which  is 
also  solved  numerically  to  obtain  the  flow  distribution.  A 
tridiagonal  matrix  algorithm  (TDMA)  was  implemented  to 
obtain  pore  velocity  profiles  and  the  electroosmotic  drag  co¬ 
efficient  in  the  capillary  pores. 


4.  Results  and  discussions 

Non-dimensional  potential  and  Navier-Stokes  equations 
were  solved  simultaneously  to  obtain  the  electroosmotic  flow 
and  to  identify  its  flow  characteristics  in  a  PEM  pore  structure 
under  a  variety  of  design  and  operating  conditions.  All  of  the 
simulations  are  based  on  fully  hydrated  Nafion  membranes 
under  the  condition  that  the  water  content  of  the  membrane  on 
the  anode  side  is  sufficiently  high  to  replenish  the  water  and 
maintain  the  membrane  hydration.  The  membrane  thickness 
is  considered  to  be  inflated  by  30%  to  account  for  the  hydra¬ 
tion.  The  parameters  studied  and  the  cell  physical  properties 
are  summarized  in  Table  1 . 


Table  1 

Parameters  and  properties  used  in  the  present  PEM  simulations 


Quantity 

Value 

Cell  temperature 

80  °C 

Pore  water  viscosity 

3.6  x  10_4kgm_1  s-1 

Pore  water  density 

980  kg  m-3 

Cell  voltage 

0.70  V 

Fixed  charge  site  concentration  [31] 

1.2  x  10-3  mol  cm-3 

Wet  membrane  thicknessa 

Nafion  117 

2.31  x  10-4  m 

Nafion  115 

1.65  x  10-4  m 

Nafion  112 

6.60  x  10“5  m 

Inflated  thickness. 
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4.1.  Estimation  off  potential 

In  order  to  solve  the  governing  equations,  one  needs  to 
know  the  wall  potential  (f  potential).  To  the  best  knowledge 
of  the  authors,  this  quantity  has  not  been  yet  reported  for 
the  membrane  in  PEM  fuel  cell  under  operating  conditions. 
Therefore,  in  the  first  set  of  calculations,  the  f  potential  was 
varied  over  a  reasonable  range  (from  —400  to  — 100  mV)  and 
the  resulting  electroosmotic  drag  coefficients  were  compared 
with  the  already  reported  values  in  the  literature. 

Fig.  4a  depicts  the  variation  of  the  non-dimensional  charge 
density,  p*,  as  a  function  of  the  non-dimensional  pore  radius, 
r*,  for  three  different  f  potential  values  and  zero  pressure 
difference  between  the  anode  and  cathode  ends.  As  expected, 
the  sulphonic  wall  attracts  free  protons  creating  a  large  charge 
density  at  the  boundary.  The  charge  density  decreases  rapidly 
away  from  the  wall  and  levels  off  at  the  pore  center.  The 
resulting  non-dimensional  potential  and  velocity  distribution 
are  shown  in  Fig.  4b  and  c,  respectively.  As  seen  in  Fig.  4b,  the 
electric  potential  has  a  minimum  at  the  wall  r*  =  1  (i.e.  the  f 
potential),  and  increases  rapidly  towards  the  pore  center.  The 
effect  of  the  cell  electric  field  is  thus  maximum  in  the  near 
wall  region  where  p*  is  at  the  highest.  From  Fig.  4c,  it  is  seen 
that  a  non-zero  velocity  (i/g)  exists  at  the  wall  which  increases 
slowly  in  the  direction  of  the  pore  centerline  such  that  a  plug- 
type  flow  distribution  prevails  in  the  pore.  The  larger  the  f 
potential  (absolute  value),  the  stronger  the  driving  forces  (by 
increasing  p*),  thus  the  larger  the  electroosmotic  flow  in  the 
pore. 

The  volumetric  flow  rate  arising  from  the  electroosmotic 
flow  in  the  capillary  pore  was  calculated  by  integrating  the 
velocity  distribution  across  the  pore  area.  Fig.  5a  displays 
the  variation  of  the  flow  rate  as  a  function  of  the  f  potential. 
The  electroosmotic  drag  coefficient  can  then  be  calculated 
from: 


K  =  NU20  =  (P_\  NAlfumirydr 

Dra§  iVH+  VM/ Water  f  NH+(r)upo(r)r  dr 

where  Na  is  the  Avogadro  number,  Mis  the  molecular  weight, 
VH+  is  the  number  density  for  protons,  and  ueo  is  the  velocity 
distribution  due  to  the  electroosmotic  flow. 

The  variation  of  Kj)rSLg  as  a  function  of  f  potential  is  shown 
in  Fig.  5b.  As  seen  from  this  figure,  the  most  reasonable 
f  potential  is  the  one  giving  rise  to  Kj)mg  =  1,  comparable 
with  the  previous  studies  summarized  in  Section  1 .  Therefore, 
the  corresponding  value  of  f  =  —0.327  V  is  taken  for  the 
current  parametric  study  calculations  which  are  presented  in 
the  following  sections. 


(a)  r* 


(b)  r 


10 
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(C) 


Fig.  4.  Effect  of  zeta  potential  on  the  distribution  of  (a)  the  charge  den¬ 
sity,  (b)  the  potential,  and  (c)  the  resulting  electroosmotic  velocity  through 
the  fully  hydrated  Nation  117  membrane  with  the  parameter  values:  if/fr  = 
(-13.15,  —6.57,  —3.29),  P*  =0,M*  =  (-3.48,-1.74,  —0.87),  and  E*  = 
(-1.75,  -3.5,  -7.0). 


4.2.  Effect  of  pore  diameter 

There  exists  a  pore  size  distribution  in  the  Nafion  mem¬ 
brane.  Therefore,  it  would  be  useful  to  investigate  the  effect  of 
pore  diameter  on  the  electroosmotic  flow.  Theoretically,  pore 
diameter  not  only  affects  the  diffuse  layer  thickness  {ic\a  in 


Eq.  (5)),  but  also  influences  the  strength  of  the  viscous  forces 
(M*)  in  Eq.  (14).  Fig.  6a  shows  that  the  larger  the  pore  di¬ 
ameter,  the  lesser  the  charge  density  becomes  towards  the 
pore  centerline.  Also,  the  larger  source  term  in  Eq.  (5)  re¬ 
sults  in  a  larger  potential  at  the  pore  centerline,  as  shown  in 
Fig.  6b.  The  larger  diameter  further  enhances  electroosmot- 
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Fig.  5.  The  variation  of  (a)  a  single  pore  flow  throughput  and  (b)  the  elec- 
troosmotic  drag  coefficient  as  a  function  of  the  £  potential  for  the  fully 
hydrated  Nation  117  membrane. 

ically  induced  flow  velocity  in  the  pore  structure,  evident  in 
Fig.  6c. 

With  no  pressure  difference,  Fig.  7a  shows  how  the  elec- 
troosmotic  flow  through  a  single  pore  increases  non-linearly 
with  the  pore  size.  Variation  of  the  /^Drag  with  the  pore  size 
given  in  Fig.  7b  reveals  that  this  parameter  is  almost  linearly 
changing  with  the  pore  diameter. 

4.3.  Effect  of  electric  field  strength 

Both  membrane  thickness  and  the  cell  voltage  affect 
the  electric  field  strength  across  the  electrolyte  membrane. 
Although  the  potential  distribution  is  unaffected,  the  elec- 
troosmotic  velocity  increases  linearly  with  E  (Eq.  (14)). 
This  linear  dependency  has  been  experimentally  observed 
for  electric  field  strengths  of  up  to  3  x  107  V  m-1  [30]. 

If  the  cell  voltage  remains  constant,  will  inversely  de¬ 
crease  with  the  membrane  thickness  resulting  in  a  signif¬ 
icantly  lower  electroosmotic  flow  (or  higher  water  content) 
for  thicker  membranes  such  as  Nation  1 17  as  shown  in  Fig.  8. 


Fig.  6.  Effect  of  pore  diameter  on  the  distribution  of  (a)  the  charge  density, 
(b)  the  potential,  and  (c)  the  resulting  electroosmotic  velocity  for  the  fully 
hydrated  Nation  117  membrane  with  the  parameter  values:  =  —10.75, 

P*  =  0,  M*  =  (-0.46,  -2.84,  -7.28),  and  E*  =  -2.14. 

However,  this  could  be  balanced  with  the  lower  back  diffusion 
and  higher  ohmic  resistance  usually  involved  with  thicker 
membranes.  On  the  other  hand,  is  unaltered  as  the  po¬ 
tential  distribution  is  unaffected  by  the  membrane  thickness. 
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Fig.  7.  The  variation  of  (a)  a  single  pore  flow  throughput  and  (b)  the  elec- 
troosmotic  drag  coefficient  as  a  function  of  the  pore  diameter  for  the  fully 
hydrated  Nation  117  membrane. 

4.4.  Effect  of  pressure  differential 


Fig.  8.  The  variation  of  electroosmotic  flow  as  a  function  of  the  PEM  thick¬ 
ness,  L,  for  the  fully  hydrated  Nation  membrane  with  the  parameter  values: 

=  —10.75;  P*  =  0;  M*  =  -9.96  for  Nation  112, -3.98  for  Nation  115, 
and  —2.84  for  Nation  117. 

With  the  electroosmotic  flow  rates  determined,  the  APs  cal¬ 
culated  from  Eq.  (19)  are  extremely  large  and  impractical 
to  apply.  This  is  due  to  the  very  small  pore  diameters  in  the 
membrane.  Fig.  9  indicates  that  a  pressure  difference  of  about 
19,000  atm  is  required  to  stop  completely  the  electroosmotic 
flow  from  anode  to  cathode  in  Nation  117  membrane.  How¬ 
ever,  in  practice,  application  of  a  few  atm  pressure  differences 
has  provided  a  better  cell  performance  due  to  reduced  drying 
on  the  anode  side  of  the  membrane. 

4.5.  Effect  of  sulphuric  acid 

Preparation  of  polymer  electrolyte  membrane  is  a  key  fac¬ 
tor  for  the  fuel  cell  performance.  During  the  preparation,  the 
solid  electrolyte  integrated  with  the  catalyst  layers  and  elec- 


Net  water  transport  in  the  membrane  arises  from  three 
mechanisms:  electroosmotic  drag,  convection,  and  diffusion. 
The  effect  of  the  first  two  mechanisms  on  the  membrane  con¬ 
ductivity  is  shown  in  Eq.  (12).  The  momentum  equation  is 
linear  and  one  can  decompose  the  velocity  field  into  two  parts: 

=  ul0  +  u\P  (18) 

where  WEO  and  u*AP  correspond  to  the  non-dimensional  elec¬ 
troosmotic  and  pressure-driven  capillary  flow  velocities,  re¬ 
spectively. 

If  the  electroosmotic  flow  is  to  be  stopped  exclusively  by 
applying  a  pressure  difference  across  the  membrane,  one  may 
use  the  Poiseuille  equation  for  fully  developed  flow  in  circular 
tubes  to  estimate  the  required  pressure  as  follows: 


A  P  = 


SjlLQ 
na 4 


(19) 


Fig.  9.  Counterbalance  of  electroosmotic  and  pressure-driven  flows  in  a  fully 
hydrated,  fully  developed  Nafion  membrane  pore  with  the  parameter  values: 
^  =  -10.75,  P*  =  -144,  M*  =  -2.84,  and  E*  =  -2.14. 
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trodes  is  cleaned  by  immersing  in  boiling  hydrogen  peroxide 
solution.  This  is  followed  by  boiling  it  in  sulphuric  acid  to  en¬ 
sure  a  full  protonation  of  the  sulphonate  groups.  The  process 
is  completed  by  rinsing  the  membrane  in  boiling  de-ionized 
water  to  remove  any  remaining  acid.  It  is  therefore  possible 
for  some  sulphuric  acid  traces  to  escape  the  rinsing  process. 
In  some  applications,  sulphuric  acid  is  intentionally  added 
to  the  membrane  to  enhance  its  conductivity  properties.  In 
this  section,  the  effect  of  the  presence  of  sulphuric  acid  in  the 
membrane  on  the  electroosmotic  flow  is  investigated. 

Ionization  of  sulphuric  acid  in  the  membrane  modifies 
the  Poisson-Boltzmann  equation,  Eq.  (5),  by  influencing 
the  charge  density  distribution.  With  extra  H+  and  SCU2- 
ions  present,  the  non-dimensional  net  charge  density,  p*,  be¬ 
comes: 

p*  =  exp(— r/r*)  +  2n2,°°  [exp(— x/r*)  —  exp(2^*)l  (20) 

n0 o 

where  ^2,00  is  the  sulphuric  acid  bulk  number  concentration. 
The  contributions  from  the  acid,  H+  and  SO42-  ions,  are  rep¬ 
resented  by  the  second  and  third  terms  in  Eq.  (20),  respec¬ 
tively.  The  modified  Poisson-Boltzmann  equation  becomes: 

1 A  (r*  A!) 

r*  dr*  \  dr*  ) 

=  -^(/ci<2)2exp(-^*)  +  (/ K2a )2  [exp(2^*)  exp(-^*)] 

(21) 

where  K2  is  the  same  as  k\  except  nQ 0  replaced  by  722,00-  The 
change  in  the  potential  distribution  due  to  the  presence  of  sul¬ 
phuric  acid  also  changes  the  source  term  in  the  momentum 
equation  (e.g.  in  Eq.  (12))  and  modifies  the  electroosmotic 
flow  distribution  in  the  membrane.  Due  to  the  lack  of  exper¬ 
imental  or  theoretical  data,  it  is  assumed  that  acid  presence 
does  not  influence  the  f  potential  appreciably. 

The  effect  of  sulphuric  acid  concentration  on  the  elec¬ 
troosmotic  flow  was  studied  by  varying  the  acid  concentra¬ 
tion  from  0  to  2M.  Fig.  10a  shows  that  the  electroosmotic 
flow  increases  steadily  with  acid  concentration.  The  water 
transport  due  to  electroosmosis  is  almost  doubled  at  2  M  acid 
concentration  compared  with  that  of  non-doped  membrane. 
The  value  of  T^Drag  however  falls  steadily  with  acid  concen¬ 
tration,  as  shown  in  Fig.  10b.  This  is  because  a  larger  number 
of  protons  is  now  taken  into  account  for  the  calculations  in 
Eq.  (17). 

The  variation  of  the  electroosmotic  flow  with  acid  concen¬ 
trations  can  be  attributed  to  the  changes  in  ion  dispersions  and 
the  resulting  net  charge  density  distribution  in  the  membrane. 
Fig.  11  shows  how  the  H+  (from  both  PEM  and  acid)  and 
S042-  ions  are  distributed  in  the  pore  at  1  M  acid  concentra¬ 
tion.  As  expected,  SO42-  ions,  now  present  in  the  membrane, 
are  repelled  from  the  negative  wall  and  at  the  same  time,  more 
protons  are  clustered  in  the  close  vicinity  of  the  wall.  As  a 
result,  the  net  charge  distribution  is  affected  which  in  turn 
modifies  both  the  potential  and  the  electroosmotic  flow  dis- 


Fig.  10.  (a)  Electroosmotic  flow  and  (b)  electroosmotic  drag  coefficient  as  a 
function  of  sulphuric  acid  molarity  in  fully  hydrated  Nation  117  membrane. 


Fig.  1 1 .  Number  density  distributions  for  PEM  protons  and  the  contributions 
from  sulphuric  acid  (1  M)  as  a  function  of  the  pore  radius  for  fully  hydrated 
Nation  117  membrane  with  the  parameter  values:  =  —10.74,  P*  =  0, 

M*  =  -2.84,  and  E*  =  -2.14. 
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tributions  as  shown  in  Fig.  12.  The  extra  protons  near  the 
wall  increase  the  charge  density  there,  and  additional  SO42- 
in  the  centerline  reduces  the  charge  density  in  central  regions 
as  indicated  in  Fig.  12a.  This  causes  the  potential  to  increase 


Fig.  12.  Effect  of  sulphuric  acid  concentration  on  the  distribution  of  (a)  the 
charge  density,  (b)  the  potential,  and  (c)  the  resulting  electroosmotic  velocity 
through  the  fully  hydrated  Nation  117  membrane  with  the  parameter  values: 
x[r^  =  -10.74,  P*  =  0,  M*  =  -2.84,  and  E*  =  -2.14. 


sharply  in  the  near  wall  region  and  drop  in  the  center.  Nu¬ 
merical  results  show  that  the  wall  effects  are  stronger  at  all 
acid  concentrations  as  displayed  in  Figs.  10  and  12b,c. 

5.  Conclusions 

The  Poisson-Boltzmann  and  the  Navier-Stokes  equations 
were  solved  numerically  to  determine  the  electroosmotic  flow 
and  the  resulting  electroosmotic  drag  coefficient  in  an  oper¬ 
ating  PEM  fuel  cell.  Based  on  the  physical  model  and  the 
numerical  results,  following  conclusions  can  be  drawn:  (a) 
Electroosmotic  flow  through  a  single  pore  increases  non- 
linearly  with  the  pore  size.  The  resulting  ^Drag  increases 
linearly  with  the  pore  diameter,  (b)  The  membrane  thickness 
has  a  significant  effect  on  the  electroosmotic  flow.  At  con¬ 
stant  cell  voltage,  thicker  membranes  (e.g.  Nation  117)  pro¬ 
vide  more  resistance  for  the  electroosmotic  flow.  However, 
^Drag  is  not  affected  by  the  membrane  thickness,  (c)  Due 
to  the  very  small  pore  diameters  (e.g.  <50  A),  the  pressure 
difference  required  to  completely  counterbalance  the  elec¬ 
troosmotic  flow  is  very  large  and  impractical  to  apply,  (d) 
When  membrane  is  doped  with  sulphuric  acid,  numerical  re¬ 
sults  revealed  that  the  electroosmotic  flow  increases  steadily 
with  the  acid  concentration. 
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